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PAPER  NO.  21 


THE  EFFECT  OF  IN  IT  LAI  CURVATURE  ON  THE  DYNAMIC  RESPONSE 
OF  THE  SPINE  TO  AXIAL  ACCELERATION 

T.F.  Li,  S.H.  Advani ,  and  Y-C.  Lee 

Department  of  Theoretical  arid  Applied  Mechanics 
West  Virginia  University,  Morgantown,  West  Virginia 


ABSTRACT 

A  majority  of  the  studies  on  the  dynamic  response  of  the  human 
torso  have  considered  uni -axial  models  wherein  the  initial  curvature 
of  the  spine  is  ignored.  A  detailed  discrete  parameter  vertebral 
response  model  incorporating  the  variable  geometry  of  the  spine  and 
subjected  to  pilot  ejection  simulated  impact  conditions  has  been 
recently  investigated  by  Qrne  and  Liu.  In  this  work,  a  simple  con¬ 
tinuum  representation  of  the  spine  is  formulated  and  the  resulting 
boundary  value  problem  is  solved  for  the  axial  and  lateral  (bending) 
dynamic  response.  The  assumed  model  is  a  constant  cross-section, 
sinusoidally  curved,  elastic  beam  with  an  end  mass  subjected  to  an 
axial  acceleration  at  the  other  end.  The  effects  of  transverse  shear 
and  rotational  inertia  are  ignored  in  the  model.  The  equation  govern¬ 
ing  axial  displacement  is  a  non-homogeneous  wave  equation  subjected 
to  non-homogeneous  boundary  conditions.  The  governing  approximate 
equation  for  the  lateral  deflection  is  a  non-linear  second  order  dif¬ 
ferential  equation  with  variable  coefficients .  Short  time  solutions 
for  these  equations  are  obtained  to  demonstrate  the  effect  of  initial 
curvature  on  the  spinal  dynamic  response.  Numerical  results  indicate 
that  the  dynamic  bending  stress  is  significant  in  comparison  to  the 
axial  dynamic  stress. 
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LIST  OF  SYMBOLS 


Dimensional  Quantity 


Physical  interpretation 


A 

c  =  ^S77 

E 

I  -  Ar2 

M 

r 


effective  spinal  cross-sectional  area 

compressional  wave  speed 

instantaneous  Young's  modulus 

principal  moment  of  inertia 

concentrated  head  and  upper  torso  mass 

effective  radius  of  gyration  about 
spinal  bending  axis 

lumped  effective  torso  and  spine  mass 
density 


Non-dimensional  To  convert  to  Physical  Interpretation 

Quantity  dimensional  form 
(~)  multiply  by 


a 

c2/: 

L 

r 

P 

AE 

t 

r/c 

u 

r 

X 

r 

y(x,t) 

r 

Y(t) 

r 

y0(x) 

r 

Y0 

X  =  M/uAr  =  X  L 

r 

forcing  acceleration 
effective  spinal  length 
axial  force 
time 

axial  column  displacement 
axial  co-ordinate 
total  column  bending  displacement 
time  function 

initial  column  bending  displacement 
maximum  column  eccentricity 
mass  parameter 
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INTRODUCTION 


Crew  member  protection  from  hostile  aerospace  environments  is  a 
biomedical  engineering  problem  of  grave  concern.  Rheological  and 
structural  response  models  of  human  body  system  components  have  re¬ 
ceived  particular  attention  towards  defining  limiting  injury  thresh¬ 
olds  associated  with  the  governing  mechanism/ s)  of  injury.  Of  special 
interest  is  the  modelling  of  the  vertebral  column  response  to  transient 
headward  accelerations  along  the  spinal  axis  (+  G  ejection  mode). 
Reported  data  on  vertebral  fractures  resulting  from  pilot  ejection 
reveals  that  a  majority  of  these  fractures  occur  between  T8  and  LI. 
Several  discrete  parameter  and  continuum  models  of  the  human  torso 
ranging  in  complexity  and  scope  have  established  the  desirability  of 
analytical  representation  of  the  response  variable  defining  injury. 

A  review  of  pertinent  investigations  can  be  found  in  studies  by  von 
Gierke  [l],  Roberts,  et  al  [2],  and  Orne  and  Liu  [3].  Selected  contri¬ 
butions  are  indicated  below. 

Uniaxial  spring  mass  characterizations  of  the  human  torso  under 
impact  have  been  examined  by  Latham  i.4] ,  Payne  [5],  Stech  [6],  and 
others  [7].  A  more  refined  model  described  by  an  eight  degree  of  free¬ 
dom  damped  spring  mass  system  has  been  studied  by  Toth  [8],  Recently, 
Orne  and  Liu  [3]  have  investigated  a  detailed  multi-mass  representation 
of  the  torso  incorporating  the  effects  of  spinal  disk  axial,  bending, 
and  shear  deformations  in  addition  to  the  variable  vertebral  geometry. 

The  discrete  parameter  models  involve  the  ( simultaneous )  solution  of 
ordinary  differential  equations )  formulated  from  the  conditions  of 
dynamic  equilibrium.  Research  on  continuum  descriptions  of  the  torso 
includes  one  dimensional  wave  propagation  models  considered  by  Hess 
and  Lombard  [9],  Liu  and  Murray  [10],  Liu  [ll],  Terry  and  Roberts  [12], 
and  Murray  and  Tayler  [13].  These  uniaxial  continuum  models  vary  in 
their  degree  of  sophistication  depending  on  the  boundary  conditions 
(head  mass),  and  constitutive  relations  (linear/non-linear ,  elastic/ 
visco-elastic).  It  is  noteworthy  that  experimental  results  and  ana¬ 
lytical  solutions  comparing  the  rectangular  pulse  response  of  an 
elastic  rod-mass  system  with  that  of  an  equivalent  spring  mass  approxi¬ 
mation  have  been  obtained  by  Seigel  and  Waser  [ll*].  Their  work  indi¬ 
cates  that  the  rod-mass  system  experiences  "significantly  larger"  forces 
for  short  pulse  duration  and/or  end  mass  magnitude. 

In  this  paper,  we  consider  a  simplified  continuum  dynamic  model 
representation  of  the  curved  spine  with  the  torso  mass  uniformly  dis¬ 
tributed  along  its  length.  The  idealized  model  is  a  constant  cross- 
section,  sinusoidally  curved,  elastic  column  with  end  mass  subjected 
to  a  uniform  acceleration  at  the  other  end.  The  influence  of  trans¬ 
verse  shear  and  rotational  inertia  is  ignored  in  the  model.  In  addition, 
the  effects  of  moments  arising  from  the  head-torso  mass  eccentricity 
and  the  external  support-restraint  system  interaction  are  not  included 
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in  the  analysis.  The  bending  stress  resulting  from  these  external  force 
and  moment  intensities  can  be  superposed  on  the  selected  basic  model 
under  study.  Experimental  work  by  Vulcan,  King  and  Nakanura  [ 15 ]  indi¬ 
cates  the  relative  importance  of  support-restraint  systems  and  head- 
torso  rotation  on  bending  stresses  in  the  vertebral  beam-column.  The 
motivation  for  the  assumed  model  stems  from  the  results  of  small  animal 
+  G  impact  experiments  conducted  with  flat  back  and  contoured  support- 
restraint  systems  [16,1?].  The  large  incidence  of  vertebral  fractures 
and  paralysis  for  the  flat  back  system  evidently  supports  the  consider¬ 
ation  of  initial  spinal  curvature. 

As  a  problem  in  theoretical  mechanics,  the  non-linear  dynamic 
response  of  a  simply  supported  column  with  sinusoidal  initial  curvature 
and  a  constant  velocity  forcing  function  at  one  end  has  been  studied  by 
Hoff  [18],  Sevin  [19] »  and  Dym  and  Rasmussen  \:0],  A  comprehensive  study 
of  the  curved  dynamic  beam  response  under  constant  velocity  end  loading 
with  combinations  of  simply  supported  and  clamped  boundary  conditions 
has  been  conducted  by  Archer  and  Das  [2l].  They  demonstrate  an  improved 
numerical  stability  of  their  finite  difference  solution  when  the  effects 
of  beam  transverse  shear  and  rotational  inertia  are  included.  Here,  the 
equations  and  associated  boundary  conditions  governing  the  axial  and 
bending  spinal  column  motion  are  formulated  and  uncoupled.  The  non- 
homogeneous  wave  equation  for  the  axial  motion  is  solved  and  an  approxi¬ 
mate  equation  for  the  bending  response  time  variable  is  obtained  by 
using  the  Ritz-Galerkin  procedure.  Short  time  solutions  for  the  spinal 
response  are  obtained  by  the  Runge-Kutta  method  to  demonstrate  the 
importance  of  initial  curvature  in  considering  the  spinal  columxi 
response . 


FORMULATION  OF  THE  BOUNDARY  VALUE  PROBLEM 


The  equations  governing  motion  of  the  basic  spinal  model  (Fig.  1) 
can  be  derived  by  use  of  Hamilton's  principle.  The  Lagrangian,  using 
this  variational  energy  formulation,  considers  the  strain  and  kinetic 
energies  of  the  column  and  the  work  done  by  the  axial  force.  The  non- 
dimensional  equations,  in  terms  of  the  coupled  generalized  co-ordinates 
y(x,t)  and  u(x,t)  with  respect  to  an  inertial ly  defined  co-ordinate 
system  are: 


(i)  An  equation  governing  the  bending  motion  of  the  column 

y""  +  (P'y'  +  Py'')  +  y  =  Yq'"  (1) 

(ii)  An  equation  governing  compressive  motion  of  the  column 


-P' 


(2) 


The  non-dimensional  axial  force  P(x,t)  in  equations  (l)  and  (2)  is 
defined  by 


p  =  -|u«  +  (y')2  -  (Yq)2]| 


(3) 
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In  the  preceding  equations  primes  and  dots  denote  differentiation  with 
respect  tc  the  non-dimensionalized  space  and  time  variables  respectiv- 
ly.  The  influence  of  column  transverse  shear  and  column,  torso-head 
rotational  inertia  is  not  considered  in  these  equations.  A  more  gen¬ 
eral  representation  of  equations  (l),  (2),  and  (3)  is  currently  under 
study  at  Technology  Incorporated  using  the  finite  difference  method  [22]. 

The  boundary  and  initial  conditions  for  the  assumed  model  are 


u(0,t)  =  jj-  at2  ,  P(L,t)  =  Au(L,t)  (ta) 

y(0,t)  =  y(L,t)  =  y"(0,t)  =  y"(L,t)  =  0  (l*b) 

u(x,0)  =  u(x,0)  =  y(x,0)  =  0  ,  and  y(x,0)=  yQ(x)  (ta) 

AXIAL  DISPLACEMENT  RESPONSE  SOLUTION 
Equations  (2)  and  (3)  can  be  re-written  in  the  form 

u  -  u"  =  f'(x,t)  (5) 

with  f(x,t)  =  |[(y')2  -  (y’)2]  . 


The  solution  to  equation  (5)  can  be  obtained  by  considering  (i) 
a  homogeneous  wave  equation  subjected  to  non-homogeneous  boundary  con¬ 
ditions  and  (ii)  a  non-homogeneous  wave  equation  with  homogeneous 
boundary  conditions.  We  therefore  write 

u  *  uH  +  Up  (6) 

where  ufi  satisfies  the  equation 

^  -  4'  - 0  (T> 

with  the  boundary  and  initial  conditions 

uH(o,t)  =  |- at2  ,  u^(L,t)  =  -AUjj(L,t )  -  f(L,t)>uH(x,0) 

=  uH(x,0)=  0 

and  Up  satisfies  the  equation 

Up  -  u£’  =  f’(x,t)  (8) 

with  homogeneous  boundary  and  initial  conditions 

Up(0,t)  =  u£(L,t)  =  Up(x,0)  =  Up(x  ,0 )  =  0  . 

We  obtain  the  solution  to  equation  (7)  in  a  manner  similar  to  that  of 
Liu  and  Murray  [10]  with  the  modification  expressed  by  the  term  f(L,t) 
in  the  boundary  conditions.  The  resulting  solution  after  use  of  the 
Laplace  transform  method  is 
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00 

Z  1 

n=l  L 


Ujj(x,t)  *  —  1  (t  -  x)2  H(t  -  x)  +  y  |  (t  -  2nL  -  x)2  H(t  -  2nL  -  x) 


-  (t  -  2nl  +  x)2  H(t-2»L.x)  ^22<'1)"l£nT°C»  ' 

n=l  m=l 


•f 

J0 


(t  -  2aL  -  x)2  H(t  -  2n L  -  x)  -  (t  -  2nL  +  x)2  H(t  -  2nL  +  x) 


,m-l  -(t-t)/X  . 
(t  -  x)  e  dt 


-  (-l)n  f (L,t  +  x  -  (2a  + 

n=1  'o  ,  . 


1)L)  -  f(L»t  -  x  -  (2n  +  1)  L  • 


.  L-i  Uiil  I  xr1  4 


where  L  denotes  the  inverse  Laplace  transform  and 

e(s)  =  (1  -  As)/(l  +  As) 

The  solution  to  the  non-homogeneous  equation  (8)  can  be  obtained 
in  the  form 


"p  "Z 

n-1,3,5 


,  nnx 

5r3in2T 


f(L»t)  sin 


f(x,t)  cos  ~  dx  sin  |£-  (t  -  t)  d t  (10) 

t 

SOLUTION  FOR  THE  BENDING  RESPONSE 

Equation  (l)  governing  the  column  lateral  response  y(x,t)  is 
simplified  by  selecting  the  initial  column  deviation  yQ(x)  from  the 
vertical  axis  to  be  sinusoidal.  Additionally,  in  view  of  boundary  con¬ 
ditions  ( Ub>)  we  assume  that  the  column  responds  in  the  first  spatial 
mode.  We  therefore  take 

y(x,0)  *  y0(x)  =  Y0  sin  (wx/L)  (11) 

y(x,t)  =  Y(t)  sin  (*x/L)  (12) 
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|  Substituting  equations  (3),  (6),  (9),  (10),  (ll)  and  (12)  in  equation 

I  (1)  and  using  the  Ritz-Galerkin  averaging  method,  we  obtain  after  con- 

j:  siderable  simplification  a  non-linear,  variable  coefficient,  second 

\  order  differential  equation  governing  Y(t)  .  It  is 

f 

\  Y(t)  + 


<r  -  y0) 


-  *o> 


Y  +  - 


|A(t) 


where 


A(t) 
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22.  Y  5^ - "2~— ■  G  (t) 

Z_.  n(n2-l6)2  n 

n=l,3 


=  0 


(13) 


,  ,2 

L 

t) 

.  2  X  .  IT 

PIT  sin^  — dx  +  — 

\  L  / 

H  L  L 

0 

P*  sin  ~  cos*^  dx 
n  L  L 


3n(t)  =  j  (Y2  -  Y2)  sin  ~  (t  -  t)  dt 


and  Pjj  is  the  axial  force  derived  from  the  homogeneous  wave  equation. 

Asymptotic  solutions  to  second  order  non-linear  differential 
equations  with  variable  coefficients  of  the  type  designated  by  equation 
(15)  have  been  investigated  by  Kuzmak  (23).  However,  these  solutions 
are  valid  for  slowly  varying  time  coefficients .  A  power  series  solution 
to  (13),  with  physical  constants  represented  by  spinal  constitutive  and 
geometric  properties  and  a  forcing  acceleration  (ff)  of  20  G,  exhibited 
numerical  instability  following  a  time  duration  of  15  milliseconds.  An 
improved  technique  of  solution  using  the  Runge-Kutta  method  in  conjunc¬ 
tion  with  an  iterative  procedure  was  finally  used  to  yield  the  transient 
model  response. 

ASSUMED  SPINAL  CONSTITUTIVE  AND  GEOMETRIC  PROPERTIES 

An  extensive  literature  reviev  revealed  that  available  data  is  in¬ 
adequate  for  characterizing  the  spinal  response  in  the  short  time  domain. 
Data  surveyed  included  results  on  spinal  compressive  wave  propagation 
experiments  and  analysis  [6,7*22],  natural  axial  frequency  data  on  spinal 
response  [1,7,24],  and  compress 'on  and  bending  tests  on  human  vertebrae 
and  disks  [22,25].  The  wide  range  of  reported  results  is  evidenced  by 
comparing  (i)  a  calculated  compressive  wave  velocity  of  100  ft/sec  ob¬ 
tained  by  Hess  and  Lombard  [9]  versus  an  experimentally  determined  spinal 
cadaveric  velocity  of  191  ft/sec  [22]  and  (ii)  a  first  resonance  axial 
frequency  of  around  10  Ha  for  the  spine-upper  torso  mass  [  l]  as  compared 
to  4U  Hz  indicated  in  another  study  (7).  The  static  compressive  proper¬ 
ties  of  vertebrae  and  disks  are  well  documented  in  the  literature.  Ex¬ 
periments  conducted  at  Technology  Incorporated  (22)  on  the  compressive 
response  of  human  vertebrae  and  disks,  using  linear  visco-elastic  theory, 
indicate  that  "the  average  initial  elastic  modulus" for  sets  of  ere  bene 
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plus  one  disk  is  7^28  psi.  The  corresponding  reported  values  for  verte¬ 
brae  and  disks  are  10,029  psi  and  2552  psi  respectively  [22].  Prelimi¬ 
nary  results  from  static  moment -curvature  tests  on  human  cadaveric 
spines  reveal  that  its  flexural  rigidity  El  ranges  from  6  x  10 3  to 
104  lb-in2 [25]. 

Based  on  the  data  and  literature  reported  in  the  preceding  para¬ 
graph,  a  compromised  set  of  constitutive  and  geometric  constants  was 
selected.  The  values  chosen  for  the  assumed  model  are: 

Effective  spinal  length  (L-U  to  Cervical  Vertebrae)  L  =  18  in 
Effective  cross-section  of  area  A  =  1.3  in2 
Effective  radius  of  gyration  about  bending  axis  r  =  0.52?  in 
Spinal  column  eccentricity  =  2  in 

Head  and  upper  torso  concentrated  mass  M  =  0.055  lb  sec2/in 
Non-dimensional  mass  parameter  A  =  M/yAL  =  0.33 
Instantaneous  elastic  modulus  E  =  10,000  psi 

The  undamped  compressive  wave  velocity  and  axial  spinal  frequency 
with  the  above  data  are  120  ft/sec  and  13.5  Hz  respectively.  The 
spinal  cross-sectional  area  and  radius  of  inertia  take  into  account 
the  added  contribution  of  the  supporting  vertebral  structure.  In 
addition,  based  on  the  "hardening"  strain  rate  characteristics  of  most 
biological  materials,  the  assumed  instantaneous  model  elastic  modulus 
is  chosen  to  be  larger  than  reported  static  values. 

DISCUSSION  AND  NUMERICAL  RESULTS 

Before  proceeding  to  illustrate  the  results  of  the  numerical  compu¬ 
tations,  a  discussion  of  the  assumed  model  and  its  inherent  limitations 
will  be  presented.  The  selection  of  a  simple  half  sine  vave  for  describ¬ 
ing  the  initial  spinal  configuration  deserves  special  mention.  Based  on 
geometric  data,  the  sine  wave  adequately  defines  the  spinal  curvature 
from  the  cervical  vertebrae  to  the  upper  lumbar  region.  The  results  of 
Orne  and  Liu  [3]  demonstrate  a  vanishing  bending  moment  for  durations  up 
to  90  milliseconds  in  the  vicinity  of  the  L-3  region,  thereby  Justifying 
the  assumed  deflection  form.  In  addition,  their  results  indicate  *.  .t 
the  axial  force  remains  relatively  constant  in  the  lumbar  region  for  a 
specified  instant.  The  continuum  model,  being  on  initial  attempt  towards 
demonstrating  the  effect  of  spinal  curvature,  neither  includes  the 
bending  moment  contribution  of  the  head-torso  mass  eccentricities  nor  the 
influence  of  the  support -restraint  system  interaction.  These  effects  can 
be  introduced  by  refining  the  model  to  a  beaa-colusn  subjected  to  exter¬ 
nal  distributed  dynamic  moments  and  lateral  forces  (incorporated  as  a 
rotational  inertia  term  ««d  external  forcing  term  q(s,t)  in 

the  beam  equation).  The  effective  mass  terms  participating  in  the  bending 
and  axial  modes  would  also  require  modification  in  this  model. 


Numerical  work  was  performed  on  an  IBM  J6r,  Model  75  computer.  The 
Eunge-Eutta  routine  was  combined  with  an  iterative  procedure  to  compute 
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the  lateral  deflection.  Figure  2  illustrates  the  old-span  lateral  de¬ 
flection  time  history  computed  from  equation  (15)  for  acceleration  load¬ 
ings  of  a  =  10  G  and  20  Q.  The  corresponding  total  compressive  axial 
displacement  is  also  indicated  in  this  figure.  This  axial  displacement 
is  a  superposition  of  solutions  obtained  from  (i)  a  homogeneous  wave 
equation  (  7  )  vith  non-homogeneous  boundary  and  initial  conditions  and 
(ii)  a  wave  equation  (8)  with  a  forcing  acceleration  term  and  homo¬ 
geneous  boundary  and  initial  conditions.  Since  the  support-restraint 
reactions  are  ignored  in  the  analysis,  the  results  are  assumed  to  be 
relevant  up  to  a  time  duration  of  ho  milliseconds.  Figure  3  compares 
the  axial  force  computed  for  the  uniaxial  and  axial-bending  response 
models  at  the  load  and  head  ends.  The  significantly  reduced  dynamic 
response  factor  at  the  forcing  end  results  from  the  curvature  terms. 

The  head  end  axial  force  is  also  reduced  and  tends  to  be  tensile  for 
larger  elapsed  times.  The  maximum  compressive  fiber  stress  time  history 
at  the  anterior  mid-span  obtained  from  the  relation  o  =  -(P/A  *  Md/I) 
is  shown  in  Fig.  U.  The  assumed  initial  deflection  form  presupposes  a 
maximum  bending  contribution  at  the  mid-span  (T-8  or  T-9)  which  increases 
with  time  in  the  interval  considered.  This  is  in  contrast  with  results 
which  indicate  a  reversal  of  sign  in  the  bending  mouent  around  this 
neighborhood  (3).  The  bending  moment  stress  contribution  due  to  the 
initial  curvature  is  about  30  percent  at  t  =  ho  milliseconds.  This 
effect  would  be  further  enhanced  if  the  moments  arising  from  the  rotation 
of  the  head  and  movement  of  the  torso  were  considered.  The  instantaneous 
elastic  response  solutions  in  the  figures  represent  upper  bounds  for  the 
visco-elastic  vertebral  model.  However,  the  early  time  response  for  the 
elastic  and  damped  models  is  almost  identical. 

A  comprehensive  discussion  of  the  mechanisms  associated  with  verte¬ 
bral  injury  has  been  presented  by  Kazarian,  et  al  (26).  Among  these  the 
anterior  lip  fracture,  the  compression  fracture,  and  the  hyperextension 
fracture  are  of  particular  interest.  Coupled  with  these  findings  is  the 
reported  constitutive  experimental  data  (27,28,29).  In  the  tests  con¬ 
ducted  by  Crocker  and  Higgins  (29),  the  intervertebral  disks  exhibited 
"hardening”  stiffness  properties  with  increasing  strain  rate.  However, 
their  maximum  compressive  velocity  rate  of  bao/sec  is  well  below  the 
corresponding  rate  encountered  in  vertebral  ejection.  Based  on  the 
general  trend  of  available  experimental  results  and  the  model  analysis 
represented  by  Figs.  2,3,  and  b  the  nechanisns  of  injury  indicated  above 
lend  themselves  to  analytical  definition.  For  example,  data  on  anterior 
lip,  compression  art  hyperextension  fractures  can  be  correlated  with 
computed  dynamic  value#  of  compressive  stress,  tensile  stress  (due  to 
excessive  bending  moment)  and/or  bending  and  axial  displacements.  The 
spinous  process  fracture  with  displacement  of  '.he  pedicle  may  be  inter¬ 
preted  by  incorporating  the  effects  of  transverse  shear  in  the  governing 
equation. 
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ATIOSS 


The  geometrically  non-linear,  continuum  model  analysed  here  answers 
Some  basic  questions  pertinent  to  the  interpretation  and  prediction  of 
vertebral  column  injury  resulting  from  dynamic  axial  loads.  Specifically, 
the  initial  spinal  curvature  introduces  a  ccup'ed  axial -lateral  response. 
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FIG.  3  Axial  Force  vs.  Time,  109 


Time  -  milliseconds 

FIG.  4  Mid-span  Maximum  Fiber  Stress 


The  axial  force  and  stress  distribution  are  significantly  influenced  by 
the  bending  motion  response  as  a  result  of  the  mechanical  energy  distri¬ 
bution  between  the  axial  and  bending  modes. 

It  is  recommended  that  spinal  disk,  constitutive  equations  valid  in 
the  impact  range  be  determined  from  high  strain  rate  compressive  tests. 
Finally,  it  is  suggested  that  more  complex  continuum  models  be  formu¬ 
lated  and  the  resulting  boundary  value  problem  numerically  studied  to 
examine  the  detailed  spinal  stress  response,  the  torso  surface  wave 
response  and  the  associated  mechanical  energies  producing  injury. 
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